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Hadronic Contribution to the muon g — 2 factor 
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The lowest order hadronic contribution to the g — 2 factor of the muon is analyzed in the framework 
of the operator product expansion at short distances, and a QCD finite energy sum rule designed 
to quench the role of the e + e~ data. This procedure reduces the discrepancy between experiment 



and theory, Aa M = a^ xp — a^ M , from Aa^ 
without increasing the uncertainty. 

PACS numbers: 

I. INTRODUCTION 



The currently accepted Standard Model (SM) prediction 
of the muon magnetic anomaly is [l[ 



28.7(8.0) x 10" 



to Aa u 



20.6(8.0) x 10" 



l M = 11659180.2(4.2)(2.6)(0.2) x 10" 



(1) 



where the first error is from the lowest order hadronic 
contribution, a^ AD,LO , the second from all other 
hadronic pieces, and the third is due to all non-hadronic 
terms. In contrast, the average experimental value is 

a i 



,EXP 



11659 208.9(6.3) x l(T iU , (2) 
leading to a discrepancy between experiment and theory 



Aa,, 



EXP 



28.7(8.0) x 10 



-10 



(3) 



which is a 3.6 a effect. Another independent analysis 
of the lowest-order hadronic contribution Q gives 
o^ M = 11 659 182.8(4.9) x 10~ 10 , which corresponds to a 
3.3 tr discrepancy between experiment and theory. 



Given that the largest uncertainty in the SM predic- 
tion comes from a^ AD,LO , and that the cross-section 
data from e + e~ — > hadrons is used to evaluate this 
contribution, it is natural to wonder whether the dis- 
crepancy Aa M could be (partly) due to some issue with 
the e + e~ data. The purpose of this paper is to ex- 
amine the impact on a^ AD,LO if the contribution from 
the e + e~ data is quenched by including theoretical in- 
put in the framework of the Operator Product Expansion 



(OPE). This procedure reduces the discrepancy Aa M 
Aa„ 



to 



— 20.6(8.0) x 10 , i.e. without increasing the er- 
ror. In the process of following this procedure it will 
be shown that there is a clear discrepancy between the 
OPE based approach and the (unquenched) e + e~ data. 
These conclusions depend on two assumptions, viz. (a) 
the validity of (global) quark-hadron duality, and (b) the 
dimension d = 2 term in the OPE is entirely due to quark 
mass insertions. The validity of these assumptions will be 
critically discussed. Their impact on the hadronic con- 
tribution to the running QED coupling, and hence the 
SM fitted Higgs boson mass, will be investigated in the 
Appendix. 



II. THE METHOD 

We begin by defining the light-quark electromagnetic cur- 
rent correlator H(s), 

IW<Z 2 ) = i J d^xe^(0\T(j^(x)j^(0))\0) 

= (gnq» - q 2 g^)H(q 2 ) , (4) 

where the electromagnetic current is j^ M (x) = 
Q fQf( x )liJ. < lf( x )^ with the sum running over the light 
quarks qj- = {u, d, s}, and Qf are the quark charges. fl(s) 
is analytic in the entire complex squared energy s-planc, 
except for a branch cut along the positive real axis. Uni- 
tarity implies the optical theorem R(s) = 127rImII(s), 
where R(s) is the normalized e + e~ cross-section. The 
lowest order hadronic contribution to the anomalous 



2 



magnetic moment of the muon is given by 

,HAD,LO _ K (s)R(s)ds, 



(5) 



where K{s) is a well-known weight function. About 92% 
of the contribution to Eq. ([5]) comes from the low-energy 
region below y/s = 1.8 GcV. Here perturbative QCD 
(PQCD) is not valid, and cross-section data from e + e~ 
experiments is normally used to perform the integration, 
followed by PQCD above y/s = 1.8 GeV Q and up to 
the charm-quark threshold. Regarding the heavy-quark 
contribution, it is even possible to calculate it entirely in 
PQCD as shown recently Q. Therefore we shall concen- 
trate on the low energy region yfs < 1.8 GcV, and discuss 
how the OPE could be used to quench the contribution 
of the e + e~ data. 

We begin by considering Cauchy's theorem in the com- 
plex s-planc 



p(s)R(s) ds 



6ni <j> 

J\s\ 



p(s)U(s)ds = , (6) 



where p(s) is an arbitrary analytic function. Next, we in- 
voke quark-hadron duality to replace n(s) by IIope(s) in 
the integral around the circle, where IIope(s) is the cor- 
relator in the framework of the OPE. Using Eq. (j6|), one 
can write the low energy piece of Eq. ([5]), a^ AD,LO (so), 
as 



~HAD,LO 



(*o) 



[K(s) ~p(s)]R(s) ds 



6iri 



\s\=s 



p(s)U pe{s) ds 



(7) 



With the low energy contribution written in this form, 
one can choose a suitable p(s) to quench the contribu- 
tion of the e + e~ data, for the price of having to perform 
a contour integral. There are a number of relevant con- 
siderations in choosing p(s). First, each power s n present 
in p(s) will involve a dimension d — 2(n — 1) term in the 
OPE. As all terms beyond dimension d = 4 arc currently 
poorly known, one should avoid them by considering only 
kernels of the form p(s) = a + bs, thus implying contribu- 
tions from only dimension d = 2, 4 condensates. Higher 
dimensional terms would only contribute through loga- 
rithmic terms at order 0(a 2 ), which are expected to be 
heavily suppressed @. The second consideration is to 
note that the data in the range 1 GeV < y/s < 1.8 GcV 
is particularly badly known. The data in this region have 
the largest uncertainties, and some of the channels have 
not even been measured yet (for example the process 
e + e _ — > KKit), and their contributions must be esti- 
mated using isospin relations. Therefore, we shall fix the 
constants a, b in p(s) by demanding that the contribution 
in the region 1 GcV < y/s < 1.8 GcV be minimized. In 
particular, that the quantity 



Max 



K(s)-p(s) 



K(s) 



(lGeV < sfs < 1.8 GeV) (8) 




FIG. 1: The magnitude of the ratio of integration kernels entering 
in Eqs. {0 and (5). The function p(s) is chosen to minimize Eq. ((SJ. 



be minimized. This uniquely fixes p(s) as 

p(s) = 4.996 x 1(T 9 - 1.432 x 10 _9 s . (9) 

One can appreciate from Fig. [T]that with this choice of 
p(s) all the data in the interval \fs = (1.0 — 1.8) GeV 
is suppressed by a factor of at least 2.5. It should be 
stressed that our final results are highly insensitive to 
the use of other fitting methods, e.g. least-square fits. 
As a final comment on the choice of p(s), and in connec- 
tion with the region y/s < 1.0 GeV, the kernel K(s) is 
approximately proportional to s~ 2 , implying that K(s) 
gets flatter with increasing s. The flatter the kernel, the 
wider the range in s allowing for a linear fit. Quenching 
data at s < 1 GeV 2 over a sufficiently wide range is not 
feasible. The other issue with attempting to suppress the 
s < 1 GcV region is that K(s) — p(s) will be larger than 
K(s) in the poorly known region y/s = 1 — 1.8 GeV, thus 
emphasizing these data. This would defeat the purpose 
of this method, so that we choose instead p(s) to satisfy 
Eq. ©. 



III. OPE 



In the framework of the OPE the correlator II(s) in Eq. 
(U]) is written as 



n( g 2 



c i 



E 

N=0 



2N+2 



(q 2 ^ 2 ) (d 2N+2 (v 2 )} , (10) 



where fi is a rcnormalization scale, and where the Wil- 
son coefficients in this expansion, C2n+2(q 2 , ^ 2 ), depend 
on the Lorentz indexes and quantum numbers of J(x) 
and of the local gauge invariant operators On built from 
the quark and gluon fields. These operators are ordered 
by increasing dimensionality and the Wilson coefficients, 
calculable in PQCD, fall off by corresponding powers of 
—a 2 . In other words, this OPE achieves a factorization 
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of short distance effects encapsulated in the Wilson coef- 
ficients, and long distance dynamics present in the vac- 
uum condensates. The unit operator / in Eq. (flO|) has 
dimension d = and Col stands for the purely PQCD 
contribution 



IIpqcd(s) 



Qt 
16tt 2 



20 

y 



log - 



0(a s 



(11) 



Here Q\ = X}/=« d s Q} = 2/3. The perturbative cor- 
relator is known up to five- loop order O[ot^) @-[l3j]. We 
use as input the latest PDG world average of the strong 
coupling g|, [lj] a s (M z ) = 0.1184 ± 0.0007, which we 
run down to \i = M T using the Mathematica package 
Rundec Since there are no gauge invariant opera- 

tors of dimension d = 2 involving the quark and gluon 
fields in QCD, it is normally assumed that nonperturba- 
tive effects start at dimension d = 4. This is supported 
by results from QCD sum rule analyses of T-lepton decay 
data 0, as well as from lattice QCD (LQCD) [l7| which 
show no evidence of d = 2 operators. There is, though, a 
PQCD term of dimension d = 2 arising from quark mass 
corrections, i.e. 



c 2 (o 2 ) 



E 



Si 

Air 2 



7 



M r 



6 + 0(a s ) 



(12) 



Although these mass corrections are negligible for the 
up- and down-quarks, this is not the case for the strange 
quark. We use the three-loop 0(m 2 /s) result from (l8j |. 
and the PDG values of the light quark masses Q, i.e. 
m u (2GeV) = (2.3 ± 0.7) MeV, m d (2GeV) = (4.8 ± 
0.7) MeV and m s (2GcV) = (95 ± 5) MeV. These val- 
ues agree within their conservative errors with the most 
recent determinations from QCD sum rules Il9fl. as well 
as from (LQCD) (2(| (for a recent review see |2l|). The 
dimension d — 4 term in (TTU)) is given by Q 



C 4 (0 4 ) = 



1 



E 

f—u,d,s 



Qr 



l 

12 



11 a s (n) 



216 



-G 2 ) 



2 a s {fi) 

3 7T 



0(a 2 s 



w/(/x)(g/g/)(//) 



(13) 



The gluon condensate is known with very large uncertain- 
ties, e.g. LQCD [H gives (^G 2 ) = 0.04(1) GcV 4 , while 
QCD sum rules using the ALEPH r-decay data yield [l6[ 
(c^G 2 ) _ o oi(i) GcV 4 . There are also determinations 
based on e + e _ cross-section data [26| affected by larger 
uncertainties. We shall use here the conservative value 



G 2 ) = 0.015(15) GeV 4 



(14) 



which is consistent with most determinations. The up- 
(down-) and strange-quark condensates are well known 
to be determined from the corrections to the Gell Mann- 
Oakes-Renner relation (22|-[24[. Using this information 
it turns out that the uncertainty in the gluon conden- 
sate by far dominates over that of the second term in 



Eq.(13). Finally, we also include the QED contribution 
to the electromagnetic correlator 



IIqed(s) 



32^2 



12- 4C3 - l0S U? 



"EM 



(15) 



In order to evaluate the integrals involving R(s) in Eq. 
(O or Eq. (|5|), one needs to collate a complete set of 
e+e~ data. We shall use the data set described in (2(|. 
The procedure is to use all available data from BABAR, 
and where these data are not available, use data from the 
most recent alternative experiment. For channels where 
there is currently no experimental data, the missing data 
was estimated using isospin relations. Systematic errors 
from the same experimental group were conservatively 
taken to be 100% correlated. It should also be noted 
that for > 2 GeV, this data set employs the BES 
inclusive data [H Hi . 



IV. RESULTS 



Using the standard kernel K(s) in Eq. ([5]), and our data 
collection up to ^/sq — 1.8 GeV, we find 



had.lo _ 



(641.7 ±6.4 Adata ) x 10 



-10 



(16) 



where the error Adata is the combined statistical and sys- 
tematic error from the e + e~ data. Evaluating a^ Ar> - LO 
with the quenched kernel, as given in Eq. ([7]) , and using 
Fixed Order Perturbation Theory (FOPT) in the contour 
integral, with yfso = 1.8 GeV, and /i = M T , we obtain 



-HAD.LO 



=(650.7 ±3.7 A data± L9a conv 
± 1.0 A a s ± 1.3 a <g»>) x HT 10 
= (650.7 ±4.5) x 10~ 10 . 



(17) 



The uncertainty Aconv is the difference between the 
O(af) and O(o^) PQCD contributions, Aa s is the er- 
ror from the strong coupling, and A(G 2 ) from the uncer- 
tainty in the gluon condensate. The total error above is 
lower than that from the standard approach in Eq. (|16D , 
primarily due to the different contribution from the data, 
which amounts to 445.2 x 10~ 10 in Eq. (fTTJ), as opposed 
to 641.7 x 10~ 10 in Eq. (fTB|). Using Contour Improved 
Perturbation Theory (CIPT) to perform the contour in- 
tegral in Eq. (O we find 



-HAD.LO 



= (649.2 ±3.7 Ad ata± L3a con 
±0.6 A o.±1.3a<g»)) x 10 



-10 



(18) 



We adopt a conservative approach treating the difference 
between both methods as another source of uncertainty. 
Considering the central value to be the average of the 
FOPT and CIPT results, with an error to be added in 
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Standard Kernel , , Modified kernel 




■ jr + 7T data 

□ non-7r + 7T~ data 



500 400 300 200 100 100 200 300 400 500 

FIG. 2: The contribution to a^ AD,LO (in units of 10~ 1C1 ) from 
different sources using the standard kernel approach Eq. (0 and 
the modified kernel approach Eq. (0. 



quadrature with the other uncertainties, we obtain 

: (650.0±3.7 A data±1.3 A CO 11 V ±0.6 A 
± 1.3 A(G 2 ) ± 0.8 A mt) x 10" 10 
= (650.0 ±4.3) x 1CT 10 , (19) 



~HAD,LO 



where Aint is the uncertainty from the integration 
method (FOPT and CIPT). This result is more precise 
than the standard one, Eq. (|16p . due to the suppression 
of the data contribution achieved by using Eq.©. The 
relative contributions to a^ AD,LO from different sources 
is shown in Fig. ([2J. 
Regarding the next-to-leading order hadronic contribu- 
tions to the muon g — 2, the same e + e~ data is used in the 
standard approach based on Eq. ([5]). For instance, using 
the kernel of || gives q had,nlo = _1 .1(6) x 10~ 10 , 
which has a different sign from the leading-order con- 
tribution. Using the same procedure as above, we find 
the following difference between results from the modified 
kernel approach, a^ AD,NLO , and the standard approach, 

„HAD,NLO 



-HAD.NLO 



HAD.NLO 



-0.2 x 10 



-10 



(20) 



where no error is given, as it is not significant given the 
uncertainty in Eq. (fT9|) . Adding this next-to-leading or- 
der result to Eq. (IT^l) and recalculating the discrepancy, 
Eq. ©, gives 

Aa M = ajf p - 5* M = 20.6(8.0) x lO" 10 , (21) 

which is a lower 2.6 a effect. 



VALIDITY OF THE METHOD 



As indicated in the Introduction, the method followed 
here to determine the leading hadronic contribution to 
the muon anomaly relies on two basic assumptions, i.e. 
(a) (global) quark-hadron duality, and (b) the absence of 



a dimension d = 2 term in the OPE, beyond well known 
PQCD quark mass corrections. We discuss their validity 
in this section. 

The OPE beyond perturbation theory in the deep Eu- 
clidean region, involving a parameterization of confine- 
ment in terms of quark and gluon vacuum condensates, is 
one of the two pillars of the QCD sum rule method. This 
allows for analytic determinations of hadronic parame- 
ters, as well as of quark masses and the strong coupling 
[30j . However, the contour integral in Eq. (JT]) involves 
an integration of the OPE expression for the correlator 
near the Minkowski axis, a potential source of difficulties 
as first pointed out in the pioneering paper [3l[. Let us 
define the difference between the exact electromagnetic 
correlator and its expression as given by the OPE as 



^DV 



n(s) - IIope(s) 



(22) 



Then the exact result for the contour integral in Eq. Q 
becomes 



l=«0 



p{s)Ii(s)ds 



p(s) [IIope(s) + A DV ] ds 



(23) 

In order to obtain Eq. (O the integral over Adv was as- 
sumed to vanish, in other words no violation of (global) 
quark-hadron duality. This is supported by a model inde- 
pendent analysis [321 ] based on the vector and axial- vector 
spectral functions from hadronic r-decays as measured by 
the ALEPH collaboration (33|. However, this is currently 
a contentious issue, as there are model dependent claims 
to the contrary. In order to determine the impact of po- 
tential duality violations on our results we consider first 
a recent quantitative model for ImADv(s) [13], which 
gives 



-Im Adv(s) 

7T 



i—u,d^s 



-K V e 
6 



' sin(ay + Pvs) 



+ -K V e- lvS sin(a' y + j3 v s) 
6 



(24) 



where the values of the various constants may be found 
in [34|]. Using this model we obtain 

6ni (p p(s) Adv ds = — 127r / p(s)Im Adv (s) ds 

J\s\=s Jso 

= -0.59(59) x lO" 10 . (25) 

The error above stems from the uncertainties in the con- 
stants Ky, a' v , ay, ^VjlPv [34f. Since no error correla- 



tions were given in [341 ] we added them in quadrature. 
This result is obviously consistent with no duality viola- 
tion. In any case, it provides evidence that a potential 
duality violation is not large enough to make Eq. p^|) 
consistent with Eq. (JT7J). 

We consider next a test involving large variations in the 
value of so, for a given integration kernel p(s). This 
kernel should remain fixed, as otherwise there would be 
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o h ~ 

2 3 4 5 6 

so (GeV 2 ) 

FIG. 3: The difference <5a M = a^ AD,LO — a^ AD,LO as a function 
of so, where a M ' is given in Eq. J7J, and a" AU ' LU in Eq. ©. 
Aa M should vanish according to Cauchy's theorem. The contour 
integral in Eq. (0 was evaluated using FOPT. The dashed (solid) 
lino is the result from using the 2010 KLOE [H (2009 BABAR 
[4p| 'l data in the e + e - — > tt + tt~ channel. 

changes to the data subsets being quenched, thus pre- 
venting any conclusion. In general, one expects the size 
of any potential duality violation to decrease with in- 
creasing so- In fact, recent high-precision BES inclusive 
data [28| in the range y^i = (2.60 — 3.65) GeV are al- 
ready consistent with PQCD. In Fig. [3] we show the 
difference Sa u = a** AD ' LO - a^ AD ' LO as a function of 
s , where a^ ADXO is given in Eq. ©, and a^ AD ' LO in 
Eq. Clearly, this difference exhibits a minimum at 

y/so = 1.8 GeV, and changes around this value would 
only increase 8a M , thus reducing the current discrepancy, 
Eq. (j3j. These results also highlight the current tension 
between the OPE and the e + e~ data, and they imply a 
problem with the OPE and/or with the e + e~ data. This 
problem can hardly be blamed on duality violations, as 
8a ^ grows with increasing so- 

As it can be appreciated from Fig. [3J the kernel Eq. © 
causes 8a ^ to rapidly increase above so ~ 5 GeV 2 , al- 
though results are still consistent within the large, and 
also increasing errors. This result for Sa^ could super- 
ficially appear in contradiction with those of Q where, 
within large errors, no discrepancy was found using ker- 
nels of the form 

Pm,n(s,So)= f 1 -^) > ( 26 ) 

where the integers m, n take a variety of values in the 
range m, n = — 1. In Table 1 we show typical re- 
sults for 8a^ t using our method, but replacing Eq. 
with the above kernel. For certain combinations of the 
powers m, n the kernel Eq. (|26[) is unable to unveil the 
underlying discrepancy. However, for the combination 
m = l,n = 0, Eq. (|26p becomes similar to our kernel, 
Eq. ((5]), thus partially uncovering the tension. Since 
different kernels emphasize different regions of the data, 
some kernels arc better suited than others to uncover this 



HI 


/) 




5a M (so) 
(10 _1 °) 








2.6 GeV 


-19(52) 


1 





2.6 GeV 


13(21) 








2.0 GeV 


25(30) 


1 





2.0 GeV 


19(11) 



TABLE I: The difference Sa^ = a M ' — a M ' using the 
integration kernel [5| Eq. 126> . instead of our kernel, Eq. PJl . 
Except possibly for the case m = 1, n = 0, this kernel is unable to 
uncover the discrepancy shown in Fig. 3. 



discrepancy. 

We now examine the second assumption, i.e. the ab- 
sence of a dimension d — 2 gauge non-invariant term in 
the OPE. There is evidence in support of this assumption 
from LQCD [13], as well as from QCD sum rule analyses 
fl6l [32l of the r-decay data from the ALEPH collabora- 
tion (33J. However, this assumption could be somewhat 
relaxed, as we only require the contribution of a d = 2 
term to be positive for our conclusions to continue to 
hold (a positive contribution to Eq. (TT9")) implies a neg- 
ative value of the d = 2 term in the OPE). There is a 
model for an effective d = 2 term proposed long ago [35[ 
in connection with a so-called phcnomcnological gluon 
mass given by 

<w»>= £ e/iiv A2 (?- 32C(3) ) . (^) 

f=u,d, s V / 

where A 2 is a tachyonic gluon mass (A 2 < 0), and £(n) 
is the Riemann zeta function. The gluon mass was 
estimated in the range [H -0.085 GeV 2 < ^A 2 < 
-0.034 GeV 2 . If this term is included in the OPE Eq. 
(fT0|) . our results Eqs. (|T7l) - (fT9|) would increase by an 
amount (3.6 — 8.9) x 10~ 10 , significantly decreasing the 
discrepancy Aet M , Eq. ([3]). 



VI. CONCLUSION 



In this paper we examined within the SM the possibil- 
ity of reducing the discrepancy between the experimental 
and the theoretical value of the muon anomaly, Eq. (|3]h 
which stands as a 3.6 a effect. We made two basic as- 
sumptions, (a) the validity of (global) quark-hadron du- 
ality, and (b) the absence of a gauge non- invariant dimen- 
sion d = 2 term in the OPE, thus achieving a quench- 
ing of the contribution of the e + e~ data in the region 
1 GeV < \fs < 1.8 GeV. Using our e + e~ data compila- 
tion, and Eq. © , we achieved a reduction of the current 
discrepancy, Eq. ©, by 8.1 x 10~ 10 , leading to Eq. ([21) . 
which is now a lower 2.6 a effect. This reduction would 
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increase further by increasing sq, or by a different choice 
of data in the two-pion channel, e.g. the 2010 KLOE data 
[38^ . In this respect, the most recent KLOE analysis [39[ 
is fully consistent with [38j . and has similar uncertainties. 
In addition, we found a clear tension between our e + e~ 
data compilation and the OPE, as shown in Fig. [3] This 
suggests a cross-section deficit in the e + e~ data. We 
conclude that neither potential duality violations nor a 
hypothetical gauge non-invariant dimension d — 2 term 
in the OPE are likely to account for this current discrep- 
ancy between the OPE and the e + e~ data. 



VII. APPENDIX 



The e + e~ cross-section data used for the determination 
of the hadronic contribution to the magnetic moment of 
the muon is also used for the calculation of the hadronic 
contribution to the running QED coupling AanAoi^Iz), 
which is given by 

AaHAD(Mz) = ^l P f-J^L- ds . (28) 

^oihad(Mz) is an important input to Standard Model 
fits to the Higgs boson mass. Currently, a value of Mh = 
91^° GeV [H is obtained fr om such fits, compared to the 
mass of a potential Higgs boson detected by the ATLAS 
Collaboration with M H = 125.3(0.6) GeV j37| or by the 
CMS Collaboration with M H = 126.0(0.6) GeV Hy. It 
has been noted [43 that modifications to the e + e~ cross- 



section data that reduce the discrepancy in the muon 
g — 2, i.e. Aa M , also increases Ao!had(Mz), and hence 
decreases the fitted Higgs boson mass. This increases 
the tension between the measured and fitted Higgs boson 
masses. It is interesting in this context to examine the 
effect of reducing the contribution of the e + e~ data to 
Aochad(Mz) in the 1 GeV < y/s < 1.8 GeV region using 
an analogue of Eq. ([7]). All that needs to be changed 

is the integration kernel K(s) — > a ^ z S ^ M \ _ s ) in Eq. 

(0. We find that the linear integration kernel p a (s) = 
0.0008490-0.0002002s quenches all the data contribution 
in the interval 1 GeV < y/s < 1.8 GeV by at least a 
factor of 6.25. The greater suppression in this case is due 
to the kernel in Eq. ([28| behaving as ~ in this region, 
rather than as ~ s~ 2 for the g — 2 kernel. It is thus less 
sloped, and hence it is better fitted by the linear function 
p a (s). Denoting by <5c<had the difference between the 
modified kernel result and the standard result we find, 
using y/s~o — 1.8 GeV, 

J+1.6xl0- 4 (FOPT) 
<5«had - | +1 3 x 1Q -4 (CIpT) (29) 

This increase in AanAD(Afz) decreases the Higgs boson 
mass from Mh — 91^3 GeV [l[ to approximately Mh ~ 
(85 - 86) GeV. 
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